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Abstract 

We consider a two - dimensional model of the pseudogap state, based on the 
scenario of strong electron scattering by fluctuations of "dielectric" (AFM, 
CDW) short - range order. We construct a system of recurrence equations 
both for one - particle Green's function and vertex part, describing elec- 
tron interaction with an external field, which take into account all Feynman 
graphs for electron scattering by short - range order fluctuations. The results 
of detailed calculations of optical conductivity are presented for different ge- 
ometries (topologies) of the Fermi surface, demonstrating both the effects of 
pseudogap formation and localization effects. These results are in qualitative 
agreement with experimental data obtained for high - temperature supercon- 
ducting cuprates. 

PACS numbers: 74.25-q, 74.25.Jb, 74.72-h, 74.20.Mn 
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I. INTRODUCTION 



One of the main problems in the physics of high - temperature copper oxide supercon- 
ductors (HTSC) remains the nature of the so called pseudogap state |1|,@], existing in the 
wide region of the phase diagram. Pseudogap anomalies in electronic properties of these 
systems are observed in a number of experiments, such as measurements of optical conduc- 
tivity, NMR, inelastic neutron scattering and angle - resolved photoemission (ARPES) 
In our opinion [Q, the preferable scenario for the pseudogap formation in copper oxides is 
based on the picture of strong scattering of current carriers on well developed fluctuations of 
"dielectric" (antiferromagnetic (AFM) or charge density wave (CDW)) short - range order, 
existing in the appropriate region of the phase diagram. This scattering is strong in the 
vicinity of a characteristic scattering vector Q = (-, -) (a - is the lattice spacing of the two 
- dimensional lattice), corresponding to the doubling of lattice period (vector of antiferro- 
magnetism) and is a "precursor" of the spectrum transformation due to the appearance of 
the long - range AFM - order. Accordingly, in this region of the phase diagram develops 
essentially non - Fermi liquid like renormalization of electronic spectrum in certain parts of 
the momentum space, close to the "hot spots" on the Fermi surface , where this surface is 
effectively "destroyed" . Direct experimental confirmation of this scenario for the pseudogap 
formation was obtained in recent ARPES - experiments on Ndi^Ceo^CuO^ ||, where the 
pseudogap renormalization of electronic spectrum was observed in the vicinity of separate 
"hot spots". 

Within this scenario of the pseudogap formation we can formulate a simplified "nearly 
exactly" solvable model, describing the main anomalies of this state ||, which takes into 
account all Feynman diagrams of perturbation theory for scattering by (Gaussian) fluctu- 
ations of short - range order with characteristic scattering vector from the vicinity of Q, 
determined by the appropriate correlation length £ @,[J. This model is based on the two - 
dimensional generalization of the pseudogap formation in one - dimension due to fluctua- 
tions of CDW - type short - range order proposed rather long ago by one of the authors ||[7| . 
Simplified variant of this two - dimensional model ( "hot patches" model) was used in Refs. 
PHTT| to describe basic anomalies of superconducting state, formed on the "background" of 
the "dielectric" pseudogap. 

In Refs. [[IH mainly one - particle properties of the model were analyzed, such as the 
spectral density and the density of states. The remarkable property of this model is the 
possibility to sum also the whole Feynman series for the two - particle problem for the vertex 



part, describing the system response to an external field (e.g. electromagnetic) P JI12| , |T3|1 ■ In 
the simplified version of the "hot patches" model such calculations for the two - dimensional 
case were performed in Ref . . The main aim of the present work is the detailed analysis 
of both the theoretical aspects of calculations of two - particle properties of the more general 
model ||,|5| and calculations of the optical conductivity for different geometries (topologies) 
of the Fermi surface, appearing for "realistic" enough spectrum of non - interacting electrons. 



II. "HOT SPOTS" MODEL. 
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1. Description of the model and "nearly exact" solution for the one - particle Green's function. 

In the model of "nearly antiferromagnetic" Fermi - liquid, which is actively used to 
explain the microscopic nature of HTSC [|l|,|16|], effective electron interaction with spin - 
fluctuations of AFM short - range order is usually described by the dynamic spin suscepti- 
bility Xq(kO °f the following form, determined by fitting to NMR data |I6| : 

v ef M^) = g 2 x,H - 1 + e(q '_Q )2 _^ (i) 

where g - is a coupling constant, £ - correlation length of spin fluctuations, Q = (ir/a, n/a) 
- vector of antiferromagnetic ordering in dielectric phase, u s f - characteristic frequency of 
spin fluctuations. Both dynamic susceptibility and effective interaction (P have a sharp 
maximum in the region of q ~ Q), leading to two "types" of quasiparticles — "hot" one 
with momenta in the vicinity of "hot spots" on the Fermi surface (Fig. P and "cold" with 
momenta around the parts of the Fermi surface around Brillouin zone diagonals [Q. This 
is due to the fact that quasiparticles from the regions around the "hot spots" are strongly 
scattered with the momentum transfer of the order of Q by spin fluctuations ([!]), while for 
quasiparticles with momenta far from the "hot spots" this interaction is small enough. 

For high enough temperatures nT ^> u> s f we can neglect the spin dynamics [|J, and use 
the static approximation in ([!]): 

^ 4 ' i + at(i)- (2) 

where A - is an effective parameter with dimension of energy. In the following, similar to 
Refs. we consider both A and £ as phenomenological parameters, to be determined 

from experiments. Considerable simplification of calculations, allowing to analyze higher 
order contributions of perturbation theory, can be achieved if we replace (||) by the model 
interaction of the following form [BJ: 

v " ,(<l) = A Y 2 + (l-a^ + te-<w 2 (3) 

where A 2 = A 2 /4. Eq. (||) is qualitatively the same as @ and quantitatively is very close 
to it in most interesting region of |q — Q| < £ , determining the scattering in the vicinity 
of "hot spots" . 

We assume the following spectrum of free ( "bare" ) quasiparticles [|J : 

£ p = — 2t(cosp x a + cospya) — At cosp x acosp y a — fi (4) 

where t - is the transfer integral between nearest neighbors, while t' - is the transfer integral 
for second nearest neighbors on the square lattice, fi - is the chemical potential. This 
expression gives rather good approximation to the results of band structure calculations 
for real HTSC - systems, e.g. for YBa 2 Cu 3 6+s we have t = 0.25eV, t' = -0.45* [§. 
Chemical potential \x is fixed by carrier concentration. In this paper we shall consider 
different characteristic values of t, t' and /i, leading to different geometries (topologies) of 
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Fermi surfaces, with no relation to any particular system, with the aim of analyzing general 
enough picture. 

Consider first - order in V e ff correction to electronic self - energy, determined by the 
simplest Feynman diagram (e n = (2n + l)irT): 

£(s TO p)=£y e// (q) 1 (5) 

q »t n Sp+q 

For large enough correlation lengths £, the main contribution to the sum over q comes from 
the region close to Q = (ir/a,7r/a). Then we can write: 

£p+ q = £ P +Q+k ~ £ p +q + v p+Q k (6) 

where v p+ q = ^g+ q - is the appropriate velocity of the quasiparticle on the Fermi surface. 
Then (|5]) is easily calculated to be: 

A 2 

ie n - £ P+Q + Hl^p+ql + \v y p+Cl \)Ksigne n 

where k = Let us stress that "linearization" of quasiparticle spectrum (|6]) in @ is done 
only for a small (of the order of vf^ 1 ), due to a large values of £, correction to an energy 
spectrum of quasiparticles close to the Fermi surface, while the form of the spectra £ p and 
£p_i_q is given by general expression (|J). 

In Ref. [0 we have performed a detailed analysis of higher - order contributions for 
S(e n p). It was shown that in case of equal signs of velocity projections (on the Fermi 
surface) and Vp + Q, as well as of and "Wp+Q? Feynman integrals for any diagram of any 
order are determined only by contributions from the poles of Lorentzians in ([|) and can be 
easily calculated [} In this case the contribution of an arbitrary diagram of the iV-th order 
over effective interaction (|3]) is given by: 

= II Tl\-- (8) 

where ^(p) = ^ P+Q and v d = |^ +Q | + |^ +Q | for odd j, while ^ (p) = ^ p and Vj = + 
for even j. Here rij - is the number interaction lines surrounding j-th Green's function in a 
given diagram and for definiteness we assumed e n > 0. 

In Ref. we have given a detailed analysis of conditions for these requirements for 
velocities in the points on the Fermi surface, connected by vector Q ( "hot spots" ) to be valid, 
and also presented some examples of appropriate geometries of the Fermi surface, which are 
realized for certain relations between transfer integrals t and t' in Eq. (HJ). Under these 
conditions expression (§) is exact, with only limitation due to the use of abovementioned 
"linearization" in Eq. @. In general case (for other relations between t and t') we use 



1 Analogous statement is valid also for a special case of velocities in "hot spots", connected by 
vector Q, being exactly perpendicular to each other H. 
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(§) as a fortunate Ansatz for the contribution of an arbitrary order, obtained by simple 
continuation over spectrum parameters t and t' to the region of interest. Even in the 
worst one - dimensional case J7|, equivalent to the case of the simple square Fermi surface, 
appearing in two - dimensions for (H) with t' = and \i = 0, the use of this Ansatz gives 



results (e.g. for the density of states) which are quantitatively very close ]T7|] to the results 
of an exact numerical simulation of this problem ||18|| . In this sense we are using the term 
"nearly exact" solution. 

In case of the use of this Ansatz Eq. (|8|) it is easily seen, that a contribution of some 
arbitrary diagram with crossing interaction lines is in fact equal to the contribution of some 
definite diagram of the same order without intersections of interaction lines 0. Thus we 
can limit ourselves only to diagrams without crossing interaction lines, taking those with 
intersections into account with the help of additional combinatorial factors, attributed to 
"initial" interaction vertices (or interaction lines) @. As a result we obtain the following 
recursion relation for the one - electron Green's function (continuous fraction representation) 
0), giving an effective algorithm for numerical computations |J: 
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i£n - £*(p) + ikv k K - £ fc+ i(£ n £ p ) 

= {Gtf(en£ P ) - £ fc+1 ( £n £ p )} -1 (9) 



S fc (s n £ p ) = A 2 - V{k) (10) 

i£ n - q*(p) + ™v k K - S fc+ i(e n 4 p ) 

In graphic form this recursion relation for the Green's function is shown in Fig. [| The 
physical Green's function of interest to us is given by G(e n £, p ) = Gk=o (s n ^ p ). In (ffl) we also 
introduced a helpful technical notation: 

Combinatorial factor: 

v(k) = k (12) 

for our case of commensurate fluctuations with Q = (ir/a, ft /a) 0, if we neglect the spin 
variables (CDW - fluctuations). If we take into account the spin structure of interaction 
within the model of "nearly antiferromagnetic" Fermi liquid (spin - fermion model of Ref. 
0]), combinatorics of diagrams becomes more complicated. In this case, spin - conserving 
scattering formally gives commensurate combinatorics, while spin - flip scattering is de- 
scribed by diagrams with incommensurate combinatorics ( "charged" random field in terms 



of Ref. Q). As a result, the recursion relation for the Green's function (10) is still valid, but 
with different combinatorial factor v(k) [ffl: 

v{k) = < % (13) 
lor even k 



Below we limit ourselves to the cases of ( |P2j) and (0), some details concerning the case of 
incommensurate CDW - fluctuations can be found in Ref. |5|-|7[]. 
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Our solution for one - particle Green's function is exact in the limit of £ — > oo, when 
it can be found also in analytic form |||]]. It is a ls° exact in a trivial limit of £ — > 0, 
when for fixed value of A interaction (||) just vanishes. For all intermediate values of £ it 
gives apparently very good interpolation, being practically exact for certain geometries of 
the Fermi surface, which are realized for definite regions of spectrum parameters in || . 

Using (H) we can easily perform numerical computations of one - electron spectral den- 
sity: 

A(Ep) = —ImG R {Ep) (14) 

7T 

which can be determined also from ARPES - experiments 0. In (O) G R (Ep) denotes the 
retarded Green's function, which is obtained by the usual analytic continuation of (|9|) from 
Matsubara's frequencies to the real axis of E. Similarly, we can perform calculations of one 
- electron density of states: 

N(E) = J2 A(Ep) = -- £ ImG R (Ep) (15) 

7T — 

P P 

Details of these calculations and discussion of results for our two - dimensional model can 
be found in Refs. [g,||. 



2. Recurrence relations for the vertex part and conductivity. 

To calculate optical conductivity we need the knowledge of the vertex part describing 
electromagnetic response of the system. This vertex can be determined using the method 



suggested for the similar one - dimensional problem in Refs. [|T^,|1^]. Arbitrary diagram 
for the vertex part can be obtained by an insertion of an external field line to the appro- 
priate diagram for the self - energy || . We have already noted that in our model we can 
limit ourselves only to diagrams with non - intersecting interaction lines with additional 
combinatorial factors v(k) in "initial" interaction vertices. It is clear then that to calculate 
vertex corrections we have to consider only diagrams of the type shown in Fig. [| Then 
we immediately obtain the system of recurrence equations for the vertex parts shown by 
diagrams of Fig. [|. To find appropriate analytic expressions consider the simplest vertex 
correction shown in Fig. |] (a). Performing explicit calculations for T = in RA - channel 
we find its contribution to be: 



ri 1)RA (ep; e + up + q)=J2 V eff (K)Gi(e^. K )G R (e + ^ p - K+q ) 

K 

= A 2 {G&(e, £x(p) + iv lK ) - G R (e + u, £i(p + q) - w x k)} 1 
A 2 G^ (e, ^(p) + iv lK )G R (e + u, &(p + q) - iv lK ) \l + 



= A*Gi(e^ p )G R (e + uj, £ p+q ) {l + ■ —p?-, , , ) (16) 



^ + 6(p)-a(p + q) 

2iv\K 

^ + ^i(p)-^(p + q) 
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where during the integral calculations we have used the following identity, valid for the free 
- electron Green's functions: 

G A M P )G R (e + ^ P+q ) = {Gi(sQ ~ G$o(e + ^ p+q )} / 77 (17) 

k ' U — <5p+q + <5p 

"Dressing" the internal electronic lines we obtain the diagram shown in Fig. (b), so that 
using the identity: 



(18) 



G a (e£ p )G r (e + c^ p+q ) = {G a (e£ p ) - G R {e + ^ p+q )} x 

1 

X u - e p+q + £ P - Ef {e + wf p+q ) + Et(eQ 
valid for exact Green's functions, we can write the contribution of this diagram as: 

J RA (ep- e + cup + q) = A 2 v(l)G A (e, £ p )G R {e + lu, £ p+q ) x 
x i 1 + T~( Z u c ( \ 2W ^R( Z f — r^lu , 1 Jf A (ep;g + ^P + q) (19) 

1 W - Cl(P + q) + fl(p) - S 2 (6 + W^p+q) + £2 ( £ £ P ) J 

Here we have assumed that interaction line in the vertex correction of Fig. |5| (b) "transforms" 
self - energies J] R ' A of internal lines into S^'^,in accordance with our main approximation 
for the self - energy (cf. Fig. 0) 

Now we can write down the similar expression for the general diagram, shown in Fig. |5| 

(c): 

J k RA (ep; e + cop + q) = A 2 v(k)G A (e, £ p )G R (e + u, £ p+q ) x 

x l 1 + — n z u f n 2 Z k R k ( z f — \mva n J ^ A (gp;g + ^p + q) (20) 

1 W - 6(P + q) + &(P) - S fc+l(^ + ^P+q) + s^ +1 Kp) J 

Then we can write the general recurrence relation for the vertex part Fig. |2| in the following 
form: 



Jj^ep; £ + ^p + q) = 1 + A 2 v(k)G A (e, £ p )G R (e + cu, £ p+q ) x 

x \ 1 + n Z ufn z f — rz^TTTri ^(ep^ + ^p + q) (21) 

w - &(p + q) + &(p) - £*r+i(e + w£ p+q ) + 



The "physical" vertex J RA (ep; e+up + q) is determined as J RA Q (ep; e+up + q). Recurrence 
procedure ( pi]) takes into account a// perturbation theory diagrams for the vertex part. For 
k — > (£ — > 00) (|2l|) reduces to the series studied in Ref. || (cf. also Ref. ||), which can 
be summed exactly in analytic form. Standard "ladder" approximation corresponds in our 
scheme to the case of combinatorial factors v(k) in (|2l|) being equal to 1 |fT~3[ . 

Conductivity of our system can be expressed [|19[] via retarded density - density response 
function x R (gu;): 



2 One of the main motivations for this trick is that it guarantees the fulfillment of an exact Ward 
identity to be discussed below. 
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a(u) = e 2 Km(-^ X R (qu) (22) 

where e - is electronic charge, 

X R (qoo) = oo {$ RA (0quo) - ® RA (OOoo)} (23) 

while two - particle Green's function Q RA (equo) is determined by the loop graph shown in 
Fig. |. 

Direct numerical computations confirm that the recurrence procedure fl2T|) satisfies an 
exact relation, which directly follows (for uo — > 0) from the Ward identity derived in Ref. 
0: 

$^(0(W) = (24) 

00 

where N(Ep) - is the density of states at the Fermi level Ep = fi. This is the main argument 
for the validity of an Ansatz used to derive Eqs. flT9p, ( GUp and (pi|). 

Finally we can write conductivity in the following symmetrized form, convenient for 
numerical computations: 

^)^s^{^^P + f)^(f, P+ f;4p-f) G 14p-f)- 

(25) 

where we have introduced an additional factor of 2 due to spin summation. 

Direct numerical computations were performed using fl2"o|), (]2"T|), (^), with recurrence 
procedure being cut for high enough k, where all E fc are were supposed to be equal to 
zero. Integration in (^) was made over the whole Brillouin zone. The "bare" electronic 
spectrum was taken from (|]). Integration momenta are made dimensionless in a natural 
way with the help of the lattice constant a, while all energies in what follows are in units 
of the transfer integral t. Conductivity is measured in units of the universal conductivity in 
two - dimensions do = %- = 2.5 10~ 4 Ohm" 1 , and the density of states — in units of 1/ta 2 . 



III. RESULTS AND DISCUSSION. 

Optical conductivity and some other characteristics of the model were computed for 
different values of parameters, determining the "bare" quasiparticle spectrum (|) and for 
fixed value of A = t. Let us consider first Fermi surfaces close to the case of half - filled 
band with \x = and t' = 0, shown (for the first quadrant of the Brillouin zone) in Fig. 
[7| (a). For = and t' — the Fermi surface is a well known simple square (complete 
nesting), so that this case is quite similar to that of one - dimensional system, considered 
in Refs. []6|, I!2| , |IB | . Our results for the real part of optical conductivity of two - dimensional 



model, for the case of spin - fermion combinatorics of diagrams and different values of 
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correlation length of AFM short - range order (parameter k = where £ is measured in 
units of lattice spacing a) are shown in Fig. ||. Qualitative form of conductivity is similar 
to that found previously for the one - dimensional model (for the case of incommensurate 
CDW - fluctuations) in Refs. [0,0. It is characterized by a sharp maximum due to 
pseudogap absorption (corresponding densities of states, demonstrating pseudogap at the 
Fermi level, are shown at the insert in Fig. ^ at u ~ 2A and also by weaker maximum 
at smaller frequencies due to Anderson localization of carriers in the random field of static 
AFM fluctuations. Localization nature of this maximum is confirmed by its transformation 
into the usual "Drude - like" peak (with maximum at uj = 0) if we perform calculations 
in the "ladder" approximation, when combinatorial factors v(k) = 1, corresponding to 
"switching off" the contribution of diagrams with intersecting interaction lines, leading to 
two - dimensional Anderson localization [0,50]. Qualitative form of conductivity in this 



case is also quite similar to that obtained previously for one - dimensional model in Ref. |L3 
The narrowing of localization peak with diminishing correlation length of fluctuations can 
be explained, as was noted in Ref. [DJ, by weakening of effective interaction ([3]) for smaller 
values of £ (for the fixed value of A), leading to general weakening of scattering rate also at 
the "cold" parts of the Fermi surface. Note that general behavior of our density of states 
and optical conductivity is in full qualitative agreement with the results of quantum Monte 
- Carlo calculations of a similar model of two - dimensional Peierls transition obtained 



recently in Ref. [21 



Consider now the case of /z = again, but "switch on" the second nearest neighbor 
transfer integral t' in @. Then we are dealing with Fermi surfaces different from the square 
one, as shown in Fig. |7| (a). At the insert in this figure we show the energy dependence of 
spectral densities ( fl4T ) in several typical points of these Fermi surfaces. It can be seen that 
these spectral densities have characteristic "non - Fermi liquid" like behavior of the type 
studied in Refs. [§]||, practically everywhere at the Fermi surface, while this surface is not 
very far from the square, despite the fact that the "hot spot" in this case is precisely at the 
crossing of the zone diagonal with the Fermi surface. Corresponding frequency dependences 
of the real part of optical conductivity are shown in Fig. [|. At the insert in this figure 
we show corresponding densities of states. It is seen that as we depart from the situation 
of complete nesting, conductivity maximum due to pseudogap transitions becomes smaller, 
while the localization peak grows and becomes sharper (in accordance with the general sum 
rule for conductivity). Note, however, that pseudogap absorption peak remains noticeable 
even in the case, when the pseudogap in the density of states is practically absent (curves 4 
in Fig. D. 

Let us now return to the case of t' = and vary chemical potential /i, obtaining Fermi 
surfaces, which are sufficiently close to square, as shown in Fig. [7| (b). Strictly speaking 
there are no "hot spots" at these Fermi surfaces at all, but the spectral density, shown at 
the insert in Fig. |7| (b) still conserves a typical pseudogap form. Corresponding dependences 



for the real part of optical conductivity are shown in Fig. [TO . 

Consider now different geometries of the Fermi surface with "hot spots", shown in Fig. 
1T| . In Figs. [12], [13] we show the real part of optical conductivity, calculated (for different 



diagram combinatorics) for two characteristic values of t' = — OAt and t' = — 0.6t for the 
chemical potential value \i = 0, when "hot spots" are on the zone diagonal (curve 5 in Fig. |ll] 
(a) and curve 4 in Fig. |TTJ (b)). We can see again that the pseudogap behavior of conductivity 
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persists even for the case, when the pseudogap in the density of states (shown at the inserts 



in Figs. |12, [13]) is practically absent. Dashed curve in Fig. [12| shows the results of the 



'ladder" approximation, demonstrating typical vanishing of two - dimensional localization 



in this approximation. The results shown in Fig. [H| demonstrate typical "smearing" of the 
pseudogap maximum of conductivity with diminishing correlation length of short - range 
order fluctuations. 

For majority of copper oxide high - temperature superconductors typical geometry of 
the Fermi surface is described by the case of t' = —OAt and fi = —1.3t [|J (as shown by 
curve 3 in Fig. [TTj (a)). Results of our calculations of optical conductivity for this case 
and for different values of the inverse correlation length k are shown in Fig. [H] (for spin - 
fermion combinatorics of diagrams). Here we also introduced an additional weak scattering 
due to inelastic processes, using the standard substitution uj — > uj + ry p2]| , which leads 
to the appearance of a narrow "Drude - like" peak in the region of uj < 7 (destruction 
of two - dimensional localization due to dephasing). It can be easily checked that for 
higher values of inelastic scattering rate 7 localization peak is "smeared" and transforms 
into the "usual" Drude - like peak in the region of small frequencies. Pseudogap absorption 
maximum becomes more pronounced with the growth of correlation length £ (diminishing 
k). In Fig. ^ we show frequency dependencies of an effective scattering rate 1/t(uj) and 
effective mass m*(uj), which can be determined from our calculations, using the generalized 
Drude formula, which is often used to fit experimental data 
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m*(uj) 



m uj An 




(26) 
(27) 



Here uj p - is plasma frequency, m - is the mass of a free electron. From Fig. 15 we can see, 

UJ 2 

that \/t(oj) (expressed in this figure in units of 4^2^) demonstrates quite typical pseudogap 
behavior for the frequency range of uj < 2 A. Note, that the density of states in this case 
possess only rather weak pseudogap || (cf . insert in Fig. [12]) . In Fig. [16] we show similar 



results for the same case (typical for HTSC - oxides), but for commensurate combinatorics 
of diagrams (CDW - type fluctuations). We see that in this case the pseudogap absorption 
maximum is almost invisible. 

From Fig. [11] (b) we can see that with the change of the chemical potential in the interval 
from fi = to = — 1.666t the Fermi surface acquires larger and larger "flat" parts, trans- 
forming for « 1.666t practically into a "cross - like". Similar Fermi surface was observed 



in ARPES - experiments on Lai^sNdo^Sro ^CuO^ [£3|,|4[]. In this case velocity projections 
in "hot spots", connected by Q = (-, -), become orthogonal. For \ijt = —1.666... topology 
of the Fermi surface change (cf. Fig. |ll] (b)) and in the whole region of fi/t < —1.666... 
these projections have the same signs, guaranteeing the exactness of our main Ansatz (0) for 
the contributions of higher order diagrams ||. It is of some interest to present the results of 
our calculations of optical conductivity also for these values of fi. These are shown (for the 
case of commensurate (CDW) diagram combinatorics) in Fig. |17], where we can follow the 
changes of localization peak in conductivity as the chemical potential passes through the 
region of topological phase transition. A weak pseudogap absorption maximum practically 
does not change at all. At the insert in Fig. [17] we show the evolution of localization peak 
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as we "switch on" inelastic scattering processes (scattering rate 7) for the case of \x = —1.81 
The transition from localization to Drude - like behavior due to inelastic dephasing is clearly 
seen. These results show, that the change in topology of the Fermi surface does not lead to 
any strong qualitative change of optical conductivity in our model. 



IV. CONCLUSION 

Our analysis shows the variety of results, which may be obtained in this model for 
different geometries and topologies of the Fermi surface, appearing for different values of 
the "bare" quasiparticle energy spectrum @. It is interesting to compare these data with 
results obtained earlier for the simplified model with "hot patches" on the Fermi surface 



n|. Pseudogap anomalies in the "hot patches" model were determined mainly by strong 
scattering on these (flat) patches only, as well as by their relative size. Accordingly, the 
localization peak of conductivity in this model was rather weak and dominating behavior at 
small frequencies was determined by Drude - like peak due to the scattering on "cold" parts 
of the Fermi surface with the scattering rate 7 (analogous to inelastic scattering rate defined 
above). Analysis of the more realistic model, given in this paper, shows that the contribution 
of localization peak may be more pronounced and in fact this peak may transform into a 
narrow enough "Drude - like" peak with the inclusion of dephasing processes. 

Probably the main deficiency of the present model is the neglect of the dynamics of fluc- 
tuations of short - range order. This approximation is justified, as was noted in Refs. 
only for high enough temperatures. However, at higher temperatures inelastic scattering, 
responsible for dephasing processes and breakdown of localization, become more important. 
The other shortcoming, as was stressed several times [0,0, is our limitation to Gaussian 
fluctuations only of short - range order, which also may be justified only for high enough 
temperatures. 

Discussing the possible relation of the above results with experimental data obtained on 
real HTSC - cuprates we note, that in majority of these experiments Jl|fJ] localization peak 
was not observed at all, which may be due to the large contribution of inelastic processes 
(dephasing) at high enough temperatures used in these measurements. Peaks in optical 
conductivity at small frequencies, attributed to localization, were observed in disordered 
samples of YBaCuO in Refs. p5|j26|| . In recent experiments on NdCeCuO p7| , |28[] this peak 
was especially clearly seen. In particular, the qualitative behavior of optical conductivity 
observed in Ref. [08J for several samples of NdCeCuO with different compositions (from 
underdoped to optimally doped) is in complete agreement with our results shown above 
in Fig. |TJ, which were considered as typical for HTSC - cuprates. We conclude that the 
"hot spots" model may be successfully used for realistic enough description of anomalies of 
optical conductivity in high - temperature superconductors. 

The authors are grateful to E.Z.Kuchinskii for numerous discussions. This work was 
supported in part by the grants of the Russian Foundation of Basic Research 02-02-16031, 
CRDF No. REC-005 and the program of fundamental research of the Presidium of the 
Russian Academy of Sciences "Quantum Macrophysics" . It was also partly supported by the 
project of the Russian Ministry of Industry and Science "Studies of collective and quantum 
effects in condensed matter" . 
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FIGURES 




FIG. 1. Fermi surface with "hot spots", connected by the scattering vector of the order of 

Q = (M)- 
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FIG. 2. Diagrammatic form of the recurrence relation for the Green's function. 
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FIG. 3. General diagram for the vertex correction. 
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FIG. 6. Diagrammatic representation for two - particle response function <5 
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FIG. 7. Fermi surfaces for different values of t' and chemical potential fj,. (a) corresponds to 
H = and the following values of t'/t: - I; -0.2 - II; -0.4 - III; -0.6 - IV. (b) corresponds to t' = 
and following values of n/t: - I; -0.3 - II; -0.5 - III; -0.6 - IV. At the inserts - energy dependences 
of the spectral density for spin - fermion model for na = 0.1 at the points in the momentum space 
denoted by stars. 
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FIG. 8. Real parts of optical conductivity in spin - fermion model for the case of square Fermi 
surface (p = 0, t' = 0) for different values of inverse correlation length of short - range order rea: 
0.1 - 1; 0.2 - 2; 0.5 - 3. At the insert - appropriate densities of states. 
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FIG. 9. Real part of optical conductivity in spin - fermion model for jjl = and na = 0.1 
different Fermi surfaces, obtained from the square by "switching on" the transfer integral t'/t: 
1; -0.2 - 2; -0.4 - 3; -0.6 - 4. At the insert - appropriate densities of states. 
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FIG. 10. Real part of optical conductivity in spin - fermion model for t' = and kcl = 0.1 for 
different Fermi surfaces, obtained from the square by the shift from half - filled band. Chemical 
potentials correspond to the following values of fi/t: 0-1; -0.3 - 2; -0.5 - 3; -0.6 - 4. 
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FIG. 11. Fermi surfaces for different values of t' and chemical potential fi. (a) corresponds to 
the case of t'/t = —0.4, typical for HTSC - cuprates, and the following values of ji/t: -1.6 - 1; -1.4 
- 2; -1.3 - 3; -1.1 - 4; - 5. "Hot spots" exist for —1.6 < fi/t < 0. (b) corresponds to the case of 
t'/t = -0.6 and the following values of n/t: -1.8 - 1; -1.666 - 2; -1.63 - 3; - 4. "Hot spots" exist 
for n < 0. Dashed line - magnetic Brillouin zone. 
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FIG. 12. Real part of optical conductivity for t' /t = —0.4 and ji/t = for na = 0.1 and different 
combinatorics of diagrams: 1 - spin - fermion combinatorics; 2 - commensurate case. Dashed line 
- "ladder" approximation. At the insert - corresponding densities of states. 
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FIG. 13. Real part of optical conductivity in spin - fermion model for the case of if /t = —0.6 
and \x = for different values of inverse correlation length na: 0.1 - 1; 0.2 - 2; 0.5 - 3; 1 - 4. At 
the insert - densities of state, corresponding to 1 and 4. 
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FIG. 14. Real part of optical conductivity in spin - fermion model for t'/t = —0.4 and 
fi/t = —1.3 and different values of correlation length na: 0.05 - 1; 0.1 - 2; 0.2 - 3. Dephasing 
rate 7/t = 0.005. 
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FIG. 15. Generalized scattering rate and effective mass for the case of t'/t = —0.4 and 
fi/t = —1.3, typical for high - temperature superconductors. Parameters of the generalized Drude 
model were obtained for spin - fermion model with different values of correlation length kcl: 0.05 
- 1; 0.1 - 2; 0.2 - 3. Dephasing rate y/t = 0.005. 
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FIG. 17. Evolution of the real part of optical conductivity for commensurate case with 
t'/t = —0.6 and na = 0.2 and chemical potential changing in the vicinity of topological transi- 
tion. Different curves correspond to the following values of fi/t: -1.79 - 1; -1,77 - 2; -1.66 - 3; -1.63 
- 4. Dephasing rate 7/t = 0.01. At the insert - real part of optical conductivity for the case of 
n/t = -1.8 for different values of -y/t: 0-1; 0.01 - 2; 0.05 - 3. 
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